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Abstract 

We study a Stern-Gerlach type setup, with a quadrupole magnetic field, for neutral 
particles of arbitrary spin. The Hamiltonian is of a form proposed for joint measurements of 
incompatible observables. The measurement results are discussed, showing the limitation of 
such Hamiltonians. Some remarks are made on the relevance of covariance as a criterion for 
measurement schemes. 

1 Introduction 

The canonical form for the Hamiltonian of a measurement interaction is (operators axe careted) 

Hi = Ka' p A 0 , (1) 

where the subscripts o and p denote object and pointer, respectively. The operator A 0 is mea- 
sured, whereas a' p is conjugate to the read-out observable a p . One of the simplest measurement 
arrangements, often used as an example [1, 2, 3], is the Stern-Gerlach, where an inhomogeneous 
magnetic field effects an interaction amounting to [2, 4] 

Hi = -Kxa x , (2) 

Here the coupling constant k is proportional to the strength of the inhomogeneity. The spin 
degrees of freedom represent the object, and the spatial degrees of freedom are the “measuring 
device”. Thus we measure the spin in the x-direction, by reading out the x-momentum component. 
If we assume both p x and a x to be conserved in the absence of interactions, the Hamiltionian (2) 
causes in the Heisenberg picture (K = 1) 

Px(i) = p x {0) + na x t . (3) 

Thus, as soon as t is larger than the p x -width, the spatial part is separated into s packets, according 
to a x eigenvalue. Conversely, a read-out of p x means an accurate measurement of a x . 

In order to be able to measure two incompatible observables jointly, Arthurs & Kelly [5] 
extended the basic scheme (1) to 

Hi = K A a p A 0 + K B b' p B 0 . (4) 
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In this interaction Hamiltonian, a p and b' p are two compatible pointer system operators, with 
read-outs a p and 6 P , respectively. The operators A 0 and B 0 , not necessarily compatible, are to be 
measured. Arthurs h Kelly applied (4) to joint position-momentum measurements. They found 
that the probabilities of finding a result (a, b ) are given by 

prob(o, b) = — Tr(p 0 \a)(a\) , a = a + t6 , (5) 

7T v 7 

where |a) is a coherent or squeezed state, depending on the interaction balance ka/«b. Therefore 
the a marginal is related to the x probability distribution by a convolution 

prob(o) = f f(a — x) (x|p 0 |x) dr (6) 

J —OO 

with a Gaussian /, and analogously for b and p. This relation, providing the basis for the inter- 
pretation of b in terms of x, we have termed “non-ideality” elsewhere [6]. Sets of operators such 
as £|a)(a|, which generate probability distributions but are not orthogonal, are called positive 
operator- valued measures (POVMs) [7]. 

If we combine (4) with (2), we get a Hamiltonian like 

Hi = —K x xa x + K v ycr v . ( 7 ) 

Such a Hamiltonian may be realized using a quadrupole magnetic field, which around the origin 
in the x,y-plane satisfies 

A oc (0, 0, xy) => B (x (x, -y, 0) . (8) 

In the present paper we will investigate the properties of a measurement scheme based on the 
quadrupole Stern-Gerlach (7), where k = k x = k v . Neutral particles of arbitrary spin s will be 
used. We shall focus especially on its relation to (6) and (5). 


2 General description 

If we consider the quadrupole Hamiltonian (7), we see that its rotational symmetry immediately 
implies that J z = L z - a z is conserved. It is therefore profitable to change into a polar momentum 
representation (p, <p). Denote the eigenvalues of a z by m z . Now we can eliminate by writing the 
m z component of the state |^) as 

(m„p,(p|^) = exp[i(m ? + m*)<p] <j> m ,(p) . (9) 


The Hamiltonian, seen as an operator on the (2s + l)-dimensional spin Hilbert space, then 
becomes 


. d . . mj + a z 
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dp 


( 10 ) 


We have taken k = 2 m = 1, without loss of generality [8]. Note that OyO z is not Hermitian, but 
then neither is the overall expression (10) is Hermitian. 
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FIG. 1. Two views on the measurement results. (The Gaussian spatial part initially 
had variance A p 2 = 0.01 and L z = 0.) 

a. Typical output distribution. The density of markers indicates the momentum 
probability density at t = 1. The input state was a s = | state, with m, = — |. 

b. Matrix element (||M(p z )||) for s = \ as a function of time. The quality of the 
measurement can be characterized by the integral of the matrix element over p z > 0. 
Note that this value does not approach 1. 


We integrated eq. (10) numerically. In view of the fact that the process is intended as a spin 
measurement, we took the initial state to be a product of spatial and spin parts (denoted by 
subscripts r and a when necessary), 

I*} = \<t>)r <8 k). • (11) 

The spatial part | <f>) r was taken to be a axially symmetric Gaussian. The final state then turned 
out to be structured into a number of expanding rings, one for each |m|; only m = 0 (for integral 
spin particles) leads to a hump remaining around the origin. Fig. l.a is a typical example, where 
s = §, and we have 3 rings. In the figure the spin part was initially directed in the i-direction, 
with m x = — §, and thus we see that the distribution is peaked in the left part of the middle ring. 
More generally, states with spin initially directed in the it/-plane lead to correspondingly oriented 
distributions; only spins initially in the z-direction evolve into complete ring distributions. In the 
rings the spins are directed roughly outward, with appropriate magnitude. 

The scheme is meant to be seen as a spin measurement. Accordingly, we trace out the spatial 
variables and generate the outcome probabilities from the spin density operator p a by means of a 
POVM M(p, <p) on spin Hilbert space, 

prob- ff (p,y?) = Tr[p a M(p,<p)]-, 

M(p,<p) = r {4>\ ®i a)U{t) 1 4>) r ; U(t) =exp {-iHt) . (12) 

In the next two sections we will discuss some properties of the POVM M(p, <p). 
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3 Marginals 

The output variables of the measurement axe (p x ,p y ). Naively, one would suspect on the basis 
of (3) and (7) that a read-out of p x can be seen as a measurement of b x . Accordingly, we are 
interested in the p x maxginal of M, 

A POO A 

M(p x ) = / M(p x ,p v ) dp v . (13) 

J —CO 

Consider first the s = | case. Here the spin operators anti-commute. Thus the interaction 
Hamiltonian (7) satisfies 

SHjS 4 ' = Hr, S = /„ ® (2d x ) . (14) 

Hence, if the initial spatial state is / y -symmetric, the p x marginal of the final (p x ,p y ) distribution 
must depend only on b x . Accordingly we may write (13) as 

M(p x ) = £ f m ( Px )E m ( 0) , (15) 

where the f m axe positive functions and £ m (0) denote the projectors onto the b x eigenstates. Thus 
we get an analog of the convolution (6), reproducing the basic result (3) with additional noise terms 
that do not depend on the spin degrees of freedom [4]. As is easily verified, an analogous relation 
holds between p y and cf y . In fig. l.b we plotted the probability that a particle with m x = + 1 will 
give the measurement result p x . We see that the noise terms assure that the measurement quality 
is limited, in contrast to (3). This is necessary on account of the uncertainty principle [9, 10], as 
we are jointly measuring the two incompatible observables b x and b y . 

Nevertheless, (15) means that an unambiguous relation between the p x and b x distributions 
exists. It allows for the m x estimation from p x that we aimed at [11], albeit an imperfect one. In the 
spin- 1 case, the spin observables are Fourier-pairs [12] and therefore close analogs of the position- 
momentum pair studied by Arthurs & Kelly [5]. Accordingly, the above conclusion matches theirs. 

We might think that for s > | this result may be generalized. In fig. 2. a the diagonal elements 
of the POVM (13) in b x representation axe plotted vs. p x . We see that indeed the various m x 
values axe roughly correlated to different p x regions, as expected. Thus, it appears that from p x 
an estimate of m x can be made, just as for s = |. But there is a catch: neglecting the p 2 -term 
(strong or impulsive interaction approximation), it follows after some calculation that 

ci r\ rf 

Px(t) = Px( 0) + £*(0)[t cos 2 0 sin 2 9] 

r 

A 1 r sin rt . . A . x cosri — 1 . . 

+ <7y(0)~[ t] sin 26 + <t*( 0) sm0; (16) 

where we used the polar position representation {x y y) = (r cos 9, r sin#). Clearly p x contains extra 
spin terms that do not commute with cr x . Indeed the non-diagonal elements of M(p x ) are not zero, 
as is evidenced by fig. 2.b. Consider again (16). The desired effect, ta x , is contained in the second 
term , whereas the last two terms are the problematic ones. As the rings expand more and more, 
the i terms will decay roughly as 1/f 2 . Still, even as t — ► oo, a significant term containing a y 
remains: M and cr x are incompatible. 
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FIG. 2. The POVM M(p x ) in u x representation. (Data as in fig. 1.) 

a. Diagonal matrix elements (m|M(pj;)|m} = (— m|M(— p x )| — m). 

b. Independent non-zero non-diagonal matrix elements (m\M (—p x )\m') . 


Re(5/2,l/2) 

Im(5/2,l/2) 

Re(5/2,-3/2) 

Im(5/2,-3/2) 

Re(3/2,-l/2) 

Im(3/2,-l/2) 


But all the problematic terms in (16) are seen to contain a factor sin 26 or sin#. Denote 
reflection in the zz-plane by 7 V . The sine terms in ( p x ) vanish if we choose the spatial part of 
the initial state /^-symmetric, like we do in our calculations. But if we consider (p 2 ), terms with 
factors like <3^ sin 2 26 emerge. Such terms do not vanish, so that higher moments do not commute 
with a x . Only if s = 1, we have d 2 oc 1 so that (pi) commutes with a x after all, in accordance 
with our earlier result (15). 

But in what sense can we consider this scheme to be a measurement of cr x if s > |? Clearly 
no analog of (15) holds; only the expectation value ( p x ) is free of incompatible spin terms, so that 
we can use the measurement only to estimate the a x expectation value. Higher moments, or even 
m x probabilities, cannot be established. Thus it is a b x measurement only in the weak sense of 
expectation value estimation [11]. 


4 Covariance 


The measurement is therefore not generally a useful joint measurement. But it may have another 
use. Remember, the Hamiltonian (7) is rotationally symmetric. If we now consider the outcomes 
in polar coordinates, it is easily derived that the POVM M(p, <p) is angle covariant: 

M(p, + A<p) = R(A(p)M(p, y>)R\A<p) ; R(A<p) = exp{iA<pa z ) ■ (17) 


Covariance is a criterion that is often used to characterize classes of measurements e.g. time 
or photon phase measurements. Here we therefore speculate that M(p, <p) realizes some kind of 
spin-angle measurement. Define the projector onto the eigenstate of the operator cos 6a x + sin 9a y 
with eigenvalue m as E m (6). Then we choose as spin-angle observable [13] 



E m (6) d9 = l, 


(18) 


153 





where the summation runs over all relevant |m| values, and Cm = ^ for m = 0 and ^ otherwise. 
Accordingly, c m E m (0 ) defines a POVM that may be seen as a “spin-angle” observable [7]. We can 
now attempt to link (18) to (17). As in the previous section, this is possible for s — There a 
convolution- type relation between the realized angle measurement and the ideal (18) holds [4]. But, 
as in the previous section, this cannot be generalized to higher spins: M and 18 axe incompatible, 
although the incompatibility is generally smaller than that of fig. 2.b. Thus, although the POVM is 
angle- covaxi ant, it is not clear in what sense the “angle of spin” is measured, if at all. Analogously 
we may therefore conclude that photon-phase covariance and time covariance must also give many 
POVMs without unambiguous interpretation. Like spin-angle covariance, they axe weak criterions. 


5 Acknowledgments 

This work was sponsored by the Foundation for Philosophical Research (SWON), which is finan- 
cially supported by the Netherlands Organization for Scientific Research (NWO). 


References 

[1] R. Feynman, R. Leighton & M. Sands, The Feynman Lectures on Physics , vol. 3 (Addison- 
Wesley, Reading, Mass., 1965). 

[2] G. Ludwig, Einfiihrung in die Grundlagen der Theoretischen Physik , vol. 3 (Vieweg, Braun- 
schweig, 1976). 

[3] M. Scully, B.-G. Englert & J. Schwinger, Phys. Rev. A 40, 1775 (1989); R. Levine & R. Tucci, 
Found. Phys. 19, 175 (1990). 

[4] H. Martens & W. de Muynck, J. Phys. A 26, 2001 (1993). 

[5] E. Arthurs &; J. Kelly, Bell Syst. Techn. J. 44, 725 (1965). 

[6] H. Martens & W. de Muynck, Found. Phys. 20, 255 (1991). 

[7] C.W. Helstrom, Quantum Detection and Estimation Theory , (Academic, NY, 1976). 

[8] A suitable choice of time and position scales removes the factors in front of J^, £- and p. 

[9] H. Martens & W. de Muynck, Found. Phys. 20, 355 (1991). 

[10] P. Busch, Phys. Rev. D 33, 2253 (1986). 

[11] H. Martens & W. de Muynck, Phys. Lett. A 157, 441 (1991). 

[12] J. Schwinger, Proc. Nat. Acad. Sc. 46, 570 (1960). 

[13] The Fourier transform of a z is another candidate, but it fares even worse than (18) for the 
present (low) spin values: the incompatibility with M is much larger. 


154 



